High-visibility multi-photon interference 
of Hanbury Brown - Twiss type for classical light 
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Difference- phase (or Hanbury Brown - Twiss type) intensity interference of classical light is considered 
in higher orders in the intensity. It is shown that, while the visibility of sum-phase (NOON- type) 
interference for classical sources drops with the order of interference, the visibility of difference-phase 
interference has opposite behavior. For three-photon and four-photon interference of two coherent 
sources, the visibility can be as high as 81.8% and 94.4%, respectively. High- visibility three-photon 
and four-photon interference of space-time and polarization types has been observed in experiment, for 
both coherent and pseudo-thermal light. 
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I. INTRODUCTION 

Intensity interference, first discovered in the experi- 
ments by Hanbury Brown and Twiss (HBT) can be 
defined as the dependence of second- and higher-order 
intensity moments on certain phase delays, with the in- 
tensity itself being independent of these phase delays. 
Although the Hanbury Brown - Twiss experiment is be- 
lieved to mark the start of quantum optics, a real burst of 
interest to intensity interference started with the observa- 
tion of this effect with two-photon light [2] . The remark- 
able feature of second-order intensity interference with 
two-photon light is that its visibility can, in principle, 
be as high as 100%, while the visibility of second-order 
intensity interference with classical light cannot exceed 
50% i, 0. 

One can distinguish between two types of intensity in- 
terference, which can be called HBT-type interference 
and NOON- type interference [5]. In the HBT-type in- 
terference, the second intensity moment depends period- 
ically on the difference of the phases introduced in the 
output channels, in front of the detectors whose pho- 
tocount coincidences are registered. A typical example 
is intensity interference observed for binary stars Jp], or 
two-slit intensity interference of thermal light p3] , Q . 
It is due to this 'phase-difference' dependence that an 
HBT stellar interferometer is insensitive to atmospheric 
aberrations. NOON-type interference is observed when 
the intensity moments depend on the sum of phases . 
This type of interference has certain outstanding prop- 
erties provided that the light at the input is in an N- 
photon state: two-photon state for second-order inter- 
ference, three-photon state for third-order interference, 
and so on. The remarkable properties of NOON-type in- 
terference are super-resolution [ll| and super-sensitivity 
to phase measurements [l2|. Another important point 
is that the visibility of NOON-type interference for clas- 
sical light drops dramatically with the order of interfer- 
ence [5]. It is this fact that leads to a stronger violation 



of multi-photon Bell inequalities compared to two-photon 
ones [ill, 

However, till very recently no consideration was given 
to multi-photon HBT interference. At the same time, 
HBT interference can substitute NOON-type interfer- 
ence in important techniques called ghost diffraction 
and ghost imaging [3], [ll]. Although observed first 
with two-photon sources, ghost diffraction and imag- 
ing experiments were soon reproduced with classical 
light i, [H, 



. The only disadvantage of classical 
light with respect to two-photon interference and two- 
photon ghost imaging, compared to two-photon entan- 
gled sources, is the limited visibility, which is always be- 
low 100%. At the same time, in interference and imaging 
experiments with classical sources the visibility is inde- 
pendent of the intensity, which can therefore be arbitrar- 
ily high. In contrast, the available sources of two-photon 
light (spontaneous parametric down-conversion or spon- 
taneous four-wave mixing) should be sufficiently weak 
to provide high- visibility interference: while the visibil- 
ity is close to 100% for faint two-photon light, it drops 
with the increase in the mean number of photons per 
mode 0, 0. 

In this paper we consider higher-order intensity inter- 
ference and find, theoretically and experimentally, the 
classical limit of such interference in the third- and 
fourth-orders. The obtained values are 81.8% and 94.4%, 
respectively. We show that the visibility of HBT interfer- 
ence for classical light grows rapidly with the order of in- 
terference, in contrast to NOON-type interference. This 
property makes high-order HBT interference of classical 
sources a good candidate for ghost imaging experiments. 

The paper is organized as follows. In Section II we the- 
oretically analyze third- and fourth-order HBT-type in- 
terference of two classical light sources and derive the ex- 
pressions for the maximal visibility achievable for sources 
with coherent and thermal statistics. Section HI de- 
scribes the experiment on observing third-order inter- 
ference of two classical sources by means of coincidence 
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counting. Section IV is devoted to the polarization ana- 
logue of the observed effect. In Section V, using an alter- 
native method of processing digital images, we register 
high- visibility third- and fourth-order interference of co- 
herent and thermal light. Finally, conclusions are made 
in Section VI. 



II. THEORY 

For our consideration, we chose Young's two-slit ge- 
ometry [20]. This geometry was used in many experi- 
ments on two-photon interference and two-photon ghost 
imaging. This time, however, we consider the interfer- 
ence to be registered in the far-field zone by three de- 
tectors instead of two (Fig.l), each detector measuring 
the instantaneous intensity and triple photocount coinci- 
dences being counted by a coincidence circuit. This is the 
standard experimental technique to measure the third- 
order Glauber's intensity correlation function (ICF) [21]. 
Let A, B be point-like classical sources, having the same 
statistics, the same average intensities, and indepen- 
dently fiuctuating phases. 

If the fields of the sources are written as Ea,b = 
£^OA,Be~*^^+*'^^'^'^^^ where Eoa,b are slowly varying am- 
plitudes (for coherent beams, they are constant) and 
4'A,B{t) are the fiuctuating phases, then the instanta- 
neous intensities registered by the detectors are 

/n = /a + /b + EoAE^.se'^^-^^^^'^-^-^'^^ + c.c, (1) 

where Ia,b = \Eoa,b\'^ and (/)n = ^Au - (/>Bn, (/>5n be- 
ing the phase accumulated by the radiation of source 
S = A^B on the way to detector n = 1,2,3. Note that 
(j)n would determine the phase of the usual (first-order) 
interference pattern at point n in the absence of random 
phase fiuctuations. 

Due to the independently fiuctuating phases (j)A,B{t)^ 
there is no stable interference pattern in the far-field 
zone: time averaging makes the intensities in Eq.([T]) in- 
dependent of the phases (pn- Because of this, transverse 
displacement of the detectors (shown in Fig.l) does not 
change the intensities registered by them. However, ac- 
cording to Eq.([T]), the instantaneous intensities at points 
n = 1,2,3 are correlated or anti-correlated, depending on 
the phase difference. Correspondingly, displacement of 
the detectors will change the intensity correlation func- 
tions. Indeed, calculation of the third-order Glauber's 
correlation function for the intensities at points 1,2,3 
gives 
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- 2(cOS^i2 + COS^23 + COS^is)], 
^nm = - n, m = 1, 2, 3. 



Here, we took into account that due to time averaging, 
all terms containing the phases ^a,b(^) turn to zero. 

We see that displacement of the detectors infiuences 
only the phases (j)nm and has no effect on the intensity 
moments in ([2]), which only depend on the statistical 
properties of the sources A, B. This is because we as- 
sumed the sources to be point-like; therefore, intensity 
fiuctuations caused by each source do not depend on the 
transverse coordinates of the points 1,2,3. In Section 
HI, we will consider the case of extended sources and 
introduce the corresponding corrections into Eq.([2]). 

Passing to the normalized third-order correlation func- 
tion. 



(3) 
^123 



G 



(3) 
123 



(/l)(/2)(/3) 



, (In) = (Ia) + (Ib), 



(3) 



we obtain that 
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^123 = — + ~Y^2 ^ ^^^^^2 ^ COS^23 + COS013], (4) 



where g^'^\g^^^ are, respectively, the second-order and 
third-order normalized ICFs for each of the two sources. 
Note that not all phases are independent: (pi^, 
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The cosine sum in the square brackets of Eq. varies 
in the range from —1.5 to 3. The maximum corresponds 
to all = and the minimum, to ^12 = ^23 = 27r/3. 
To achieve these minimal and maximal values, the phases 
should vary synchronously, (j)i2 = ^23- The visibility of 
the interference pattern will then be given by the expres- 
sion 
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Because for classical light ^^^^ > [5'^^^]^ and ^^^^ > 
1 [18], with the equality holding only for coherent light, it 
follows that g'^^^ > g^'^\ and it is coherent light that pro- 
vides the largest visibility, V^H = 9/11 ^ 81.8%. Ther- 

mal radiation gives a lower visibility, V^' = 3/5 = 60%, 
which is still much higher than the corresponding value 
in the case of two-photon interference (33%). 

If the number of detectors is four, then one can measure 
fourth-order normalized correlation function, defined as 
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Writing the intensities similarly to Eq. ([T]) and calcu- 
lating fourth-order correlations, we get the expression for 
the normalized fourth-order CF in the form 
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+ cos 013 + COS 014 + COS(012 - 013) + COs(0i2 - 014) 

[^(2)12 

+ COS(013 - 014)] H [COS(012 + 013 " 014) 

+ COS(012 + 014 - 013) + COS(013 + 014 - 0i2)]. (7) 

Here, we again took into account that not all phases 
are independent. 

We note that the largest visibility is achieved in the 
case where the constant term, A = |^'^^^H-^^'^^^ + |[^'^^^]^, 
is minimal, and the amplitudes of the two oscillating 
ones, B = l^'^^^ + |[^'^^^]^ and C = |[^*^^^]^, are maximal. 
Since for classical light, g^^^ > [g^^^]^/g^^\ g^^^ > [g^'^^]^ 
and ^^^^ > 1 [18], the equalities holding only for coherent 
sources, the amplitudes A,B, C satisfy the inequalities 



A>l + l9^'\B<lg^'\C<l9^'\ (8) 

It follows that the maximal visibility corresponds to 
the case of coherent light, for which A = 1,5 = = |. 

Analysis of expression (O shows that its maximum is 
achieved at all = 0, while the minimum occurs at 
012 = f,0i3 = TT, 014 = — ^. The expression for the 
visibility becomes then 
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The visibility for coherent light is V^^l = 17/18 ^ 

94.4%. For thermal light, v}^^ = 7/9 ^ 77.8%. 

It is interesting to plot the maximal visibility values 
of this difference-phase interference, together with the 
maximal visibility values of sum-phase interference [s'], 
as functions of the order of interference. Fig. 2 shows 
these dependencies for the case of coherent light. One 
can see that while the visibility of sum-phase interference 
decreases dramatically with passing to higher orders, the 
visibility of difference-phase interference increases. 



III. EXPERIMENT WITH COINCIDENCE 
COUNTING. 

In the first series of our experiments, the third-order 
Glauber's correlation function was measured through the 
coincidence counting rate of three detectors (Fig. 3). As 
the radiation source, we used a frequency doubled Q- 
switched Nd:YAG laser with the wavelength 532 nm, 
pulse duration 5 ns, and the repetition rate 3kHz. Instead 
of two slits, which should be very precisely matched in 
width to achieve the maximum visibility of multi-photon 



interference, a single slit of width 150 /i was used, fol- 
lowed by a birefringent crystal (calcite). The crystal 
split the beam into the ordinary one and the extraor- 
dinary one; with the slit and the crystal placed between 
polarization (Glan) prisms oriented at angles ±45° to 
the plane of the crystal optic axis, this configuration was 
equivalent to two identical slits separated by a distance 
of 1.3 mm. In the far-field zone, where the interference 
pattern was formed, the radiation was attenuated using 
neutral-density filters and fed into a three-arm Hanbury 
Brown-Twiss interferometer with three photon-counting 
avalanche photodiodes and a triple-coincidence circuit 
with the coincidence resolution time 4.2 ns. Attenua- 
tion was necessary to keep the average number of photo- 
counts per pulse much less than one; otherwise, because 
of the dead-time effect, the photocount and coincidence 
rates would be measured incorrectly. Due to the gating 
of the registration electronic system, dark noise was sup- 
pressed by several orders of magnitude. In order to scan 
the interference pattern, plane-parallel glass plates with 
thicknesses 50 mm and 60 mm, respectively, followed by a 
150 /i pinholes for spatial- mode selection, were placed at 
the inputs of detectors 1 and 3. Turning the glass plates 
substituted for moving the detectors, which is usually 
performed in HBT-type experiments but is more compli- 
cated technically. By turning the plates, one could scan 
the phase of either one or two detectors within the range 
0...67r. 

In order to study multi-photon interference of sources 
with thermal statistics, we prepared pseudo-thermal light 
by means of a rotating ground-glass disc placed after 
the calcite crystal. With the disc removed, we observed 
Young's interference of two sources with coherent statis- 
tics. 

To erase first-order interference in the far-field zone we 
used an electro-optical modulator (EOM) consisting of 
four DKDP crystals. The DKDP crystals were oriented 
with their optic axes parallel to the calcite crystal axis. 
The electric field applied to the crystals changed their 
refractive indices and, as a result, the phases between or- 
thogonally polarized components of the light propagating 
through them. When an AC voltage with the frequency 
50 Hz was applied to the EOM the interference fringes 
in the far-field zone were moving. As a result, first-order 
interference in the time-averaged intensity distribution 
vanished. Such a harmonic oscillation of the relative 
phase between the ordinary and extraordinary beams led 
to the same effect as random phase fluctuations would: 
it erased the interference pattern in the time-averaged 
intensity distribution but did not influence second- and 
higher-order intensity moments. 

Due to the finite size of the slit, the two sources in 
our setup were not point-like, as it was assumed in the 
derivation of Eq.Q. In the case of coherent sources, this 
does not cause any corrections to Eq.([2|) since the normal- 
ized intensity correlation function of a coherent source is 
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equal to unity everywhere. In the case of pseudo-thermal 
sources, Eq.([2]) has to be modified if the distance between 
the detectors is comparable with or exceeds the trans- 
verse coherence length of the radiation. The transverse 
coherence length for each source, p = /A/a, is determined 
by the sizes a of the spots formed on the disc by the or- 
dinary and extraordinary beams (0.2 mm), the distance 
/ from the disc to the detectors and the radiation wave- 
length A. The first-order spatial correlation function of 
each source A.B in the far-field zone is then 



sin(27rx/p) 



2ttx/ p 
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and for thermal light, it determines all higher-order cor- 
relation functions. Then, in Eq. (|4]), the normalized ICFs 
of each of the two sources A,B g^'^\g^'^^ are not constant 
but depend on the positions ^2^1,2,3 of the detectors in the 
far-field zone: 



(2)^5^^)(x„x,)^5S = l + (7i.)^ (11) 



where 7^^ = g^^\xi — Xj), and 



1 + (712)^ + (713)^ + (723)^ + 2712713723. (12) 
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Here, the subscript 's' denotes correlation functions of a 
single source. 

Equation (pT|) is known as the Siegert relation, and 
Eq. (p!2|) is its third-order analogue. 

Then, the expression for third-order correlation func- 
tion in the case of two thermal sources becomes 
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+COS^23(723 +713712)723 

+cos(/)i3(7i3 + 712723)713], (13) 
Note that the phase differences are related to the 



coordinates 



27r6 



{xi — Xj), where b is the 



distance between the sources. 

Experimental results for the cases of coherent and 
pseudo-thermal light are shown in Fig. 4 a,b, respectively. 
Points correspond to the measured values of the normal- 
ized third-order Glauber's IGF; curves show the fit given 
by Eqs. (|4|13p . respectively, for the coherent and thermal 
cases. In accordance with the condition ^12 = — (/>32, in 
our experiment the glass plates in front of detectors 1 
and 3 were rotated synchronously, both clockwise (since 
detector 1 was in the reflected beam and detector 3 in 
the transmitted beam, this led to the opposite variation 



of the phases ^12,^32)- The obtained visibility for the 
case of coherent radiation is 74%; for the case of pseudo- 
thermal radiation, 38%. The small value of visibility in 
the case of thermal light is due to the finite transverse 
coherence length of the radiation. In each plot, we also 
show the spatial dependence of single counts for one of 
the detectors whose phase was scanned. Although the de- 
pendence is not completely fiat (the variation is caused 
by the speckle structure of the laser light and the en- 
velope of the single-slit diffraction pattern), the two-slit 
interference pattern in the intensity distribution is al- 
most completely erased. Note that the presented IGF 
was normalized to the product of the three intensities; as 
a result, the 'noisy' structure of single-photon counts did 
not infiuence the third-order interference pattern. 

From Eq. (|4|) , we can also find the interference visibility 
in the case where the third-order interference pattern is 
scanned by only one of the three detectors. This visibil- 
ity is maximal if the phase difference for the remaining 
two detectors is 7r/2. The corresponding visibility value is 
V = 1/a/2 (approximately 70.5%), which is lower than in 
the case of scanning two detectors but still considerably 
higher than in the case of the second-order interference. 
Figure 5 shows the results of g^^"^ measurement with de- 
tectors 1 and 2 fixed and detector 3 scanned. The relative 
phase of detectors 1 and 2 was aligned to be 7r/2 using 
second-order interference pattern for detectors 1 and 2. 
The observed visibility is 64%. 



IV. THIRD-ORDER POLARIZATION 
INTERFERENCE. 

The same setup, with some modification, was used to 
demonstrate the polarization analogue of the observed ef- 
fect. The modified experimental setup is shown in Fig. 6. 
By means of a quarter- wave plate placed after the EDM 
instead of the analyzer, the ordinary and extraordinary 
beams were transformed into right- and left-circularly po- 
larized beams. In the registration part of the setup, a 
polarizer was inserted in front of each detector. Because 
of the varying phase between the two beams, light in the 
far-field zone was not polarized in the first order in the 
intensity, and rotation of the polarizers did not change 
the average intensities measured by the detectors. 

At the same time, third-order correlation functions in 
the far-field zone in such a configuration do depend on the 
mutual orientations of the polarizers. Indeed, consider 
the sources A and B in Fig.l to be polarized, respec- 
tively, right- and left-circularly, and the detectors 1,2,3 
to be preceded by polarizers set at angles ^1, ^2, ^3. Then 
the instantaneous intensity at point n is given by an ex- 
pression similar to Eq.([T]), where the fields Eqa, Eqb are 
substituted by their projections onto polarization state 
selected by polarizer n, £^oA(ei?, e^), EoBi^L^^n)^ where 
the unit polarization vectors can be written in the HV 
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basis as 



= "^(^^ + ^^^)' = ^{eH - iev), (14) 

As a result, equations ([2)), (111), d) have the same form, 
with the replacement 



(15) 

If the phases are fixed, third-order and fourth-order 
intensity moments vary depending on the mutual ori- 
entations of the polarizers, although the intensities are 
constant. 

In experiment (Fig. 6), we set one of the polarizers at 
0° and rotated the other two polarizers in the opposite 
directions, with the phases (pn set at zero. According to 
Eqs.dlj), ([IS]), third-order normalized correlation function 
should vary in this case as 



dfoh = 1 + ^[cos(2^2) + cos(2^3) + cos(2^2 - 2^3)] (16) 



Its maximum should be observed at O2 = Os = 0° 
and its minimum, at O2 = —Os = 60°, and the modu- 
lation should have 81.8% visibility. Such a dependence 
was indeed obtained in experiment (Fig. 7). The visibil- 
ity achieved is 73%. Solid line shows a fit with Eq.([T6]). 
The difference between the achieved visibility and the ex- 
pected one is caused by an inaccuracy in the setting of the 
phases; for instance, dashed line shows the result of cal- 
culation with the phases set at ^12 = ^32 = 7r/6, ^13 = 0. 



V. EXPERIMENT WITH A DIGITAL 
PHOTOGRAPHIC CAMERA. 

For the measurement of fourth-order interference, in- 
stead of counting four- fold coincidences, we turned to a 
different method of measuring spatial ICFs -namely, to 
processing patterns registered by a digital camera, as sug- 
gested in Ref. [22]. The interference pattern in the far- 
field zone was recorded by a photographic camera Canon 
Powershot S2 IS. For this experiment, we used a Nd:YAG 
laser with a repetition rate 47 Hz. The exposure time 
was 1/50 s, which provided that each frame was made 
with a single laser pulse. A typical interference pattern 
recorded in one frame is shown in Fig. 8a. To accumu- 
late sufficient statistics, 500 shots were made, both for 
the coherent case and for the pseudo-thermal case. Due 
to the phase shifts introduced by the EOM, the phase of 
the interference pattern varied from frame to frame, so 
that the intensity spatial distribution averaged over all 
frames had almost no modulation (the visibility was less 



than 10%). The averaging of the intensity distribution 
and correlation functions was performed over a rectan- 
gular spatial area with the dimensions 50 pixels along y 
and 600 pixels along x (shown in Fig. 8a). This, in partic- 
ular, allowed us to eliminate the speckle structure in the 
intensity distributions (Fig. 8b), the interference modula- 
tion in which was erased to a considerable extent. The 
images were processed in the following way: first, each 
pattern like the one shown in Fig. 8a was averaged over 
50 pixels in y. This way we obtained one-dimensional 
patterns Ij{x)^j = 1,...,500. These one-dimensional 
patterns were further processed to obtain the averaged 
intensity distribution and the correlation functions. 

As it was shown above, the maximal visibility of third- 
order (fourth-order) interference patterns is achieved 
when two (three) phases 0nm are varied synchronously. 
In the triple coincidence counting measurements, this was 
done by moving two detectors simultaneously. In the 
method of digital images processing, the same result is 
achieved by measuring intensity correlations between cor- 
rectly chosen points of the image. The average intensity 
and the normalized third- and fourth-order correlation 
functions were calculated as 



n 

/(x) = (/,(x))^-^/,(x). 



{Ij{x)Ij{0)Ij{-x)) 



I{x)mi{-x) ' 
{Ij{x)Ij{0)Ij{-x)Ij{-2x)) 
I{x)I{0)I{-x)I{-2x) ' 



(17) 



where the index j numerates the frames, n is the total 
number of frames (500 in our case) and angular brackets 
denote averaging over the frames. 

Fig. 9 shows the obtained third-order (a,c) and fourth- 
order (b,d) interference patterns for coherent (a,b) and 
pseudo-thermal (c,d) sources. As expected, the distribu- 
tions in Fig 9a, c are similar to the third-order interference 
patterns registered by means of coincidence method. The 
interference visibilities achieved with the photographic 
camera are 73% and 59%, respectively. The theoreti- 
cal fit, given by Eq.(|4|) for Fig. 9a and by Eq. (p!3|) for 
Fig. 9c, is shown as a dashed line. Fourth-order interfer- 
ence fringes (Fig. 9 b,d) reveal a visibility of 93% for the 
coherent case and 81% for the pseudo- thermal case. The 
theoretical fit (using Eq.([7j)) is presented only for Fig. 9b, 
since taking into account the effects of transverse coher- 
ence in the fourth-order case leads to too bulky expres- 
sions. However, even without the fit in Fig.9d, it is clear 
that the fourth-order correlation function at the center 
is anomalously high: the theory predicts a value of 24, 
while the experiment gives a nearly twice higher value. 
The same feature, although less pronounced, is seen in 
Fig. 9c: here, the theoretical value of 6 is exceeded by 
25%. Second-order interference patterns (not presented 



6 



here) also reveal slightly excess values of ^^^^ , on the order 
of 10%. These discrepancies between theory and exper- 
iment may be explained by the non-stationarity of the 
pseudo-thermal source, caused by the intensity variation 
of the scattered light after the rotating disc. 

VI. CONCLUSION. 

In conclusion, we have considered, both theoretically 
and experimentally, higher-order HBT-type interference 
for classical light. We have shown theoretically that in 
the cases of third-order and fourth-order interference, 
the largest visibility is achieved for coherent light and 
this 'classical limit' is 81.8% in the third-order case and 
94.4% in the fourth-order one. Thermal light provides 
lower visibility of interference, which is still much higher 
than in the second-order case and reaches the values of 
60% and 77.7% for third- and fourth-order interference, 
respectively. This fact opens interesting perspectives for 
using thermal light in higher-order ghost imaging exper- 
iments [23^. Indeed, up to recently thermal light was 
supposed to be less efficient for ghost imaging than two- 
photon light, due to a lower visibility. Preliminary calcu- 
lations show that if one uses high orders in the intensity, 
this drawback of thermal light is eliminated. 

High- visibility HBT interference has been observed in 
experiment both for coherent and pseudo-thermal light. 
These results were obtained using two different meth- 
ods of correlation function measurement: by counting 
coincidences of several photodetectors and by process- 
ing a set of digital images obtained from single laser 
pulses. High-visibility third-order intensity interference 
of coherent light has been also demonstrated in polariza- 
tion measurements. The largest visibility value registered 
for coherent light was 93%, for the case of fourth-order 
interference. 

Our results demonstrate a considerable difference in 
the behavior of the HBT-type interference and NOON- 
type interference of classical light in higher orders in the 
intensity. While in the second-order intensity interfer- 
ence, classical light provides not more than 50% visibil- 
ity for both these types, the situation in higher orders 
is different. The visibility of NO ON- type interference for 
classical light is known to reduce rapidly with the order of 
interference, whereas for HBT-type interference the 'clas- 
sical visibility limit' grows with the order. Since the low 
visibility of NOON-type interference for classical light 
forms the basis for tests of two-photon and higher-order 
Bell's inequalities, it is very important to distinguish be- 
tween NOON-type and HBT-type interference in higher 
orders in the intensity. It also follows from our work 
that the mere existence of high-visibility interference in 
the third- and higher-orders in the intensity cannot be 



considered as a signature of three- or four-photon light. 
This is in contrast with two-photon interference, where 
exceeding the 50% limit is a commonly accepted criterion 
of nonclassicality. 
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FIG. 1: Young's two-slit interference experiment with the reg- 
istration of third-order intensity correlations. Radiation pro- 
duced by two sources A,B is registered by three detectors Di, 
D2, D3. Photocounts of the detectors are sent to a coinci- 
dence circuit. Transverse displacement of the detectors does 
not change their photocount rate but changes the triple coin- 
cidence counting rate. 
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FIG. 2: Maximum visibility for sum-phase (circles) \5] and 
difference- phase (squares) intensity interference as a function 
of the order of interference. 
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FIG. 3: Experimental setup. GPl and GP2, Glan prisms; 
S, single slit; EOM, electro-optic modulator; BSl and BS2, 
beam-splitters; NDF, neutral density filters; PI and P3, glass 
plates; Al-3, apertures; Dl-3, avalanche photodiodes; TCC, 
triple coincidence circuit; GGD, rotating ground-glass disc. 
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FIG. 4: Interference pattern in the normalized third-order 
intensity correlation function for (a) coherent light and (b) 
pseudo-thermal light obtained by tilting glass plates at the 
inputs of detectors 1 and 3. Empty circles show the intensity 
distribution given by the counting rate Ri of detector 1. Solid 
lines show the theoretical fit with Eq.Q (a) and Ea. (|13|) (b) 
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FIG. 5: Interference pattern in the normalized third-order 
intensity correlation function for coherent light obtained by 
tilting the glass plate at the input of detector 3. Relative 
phase between the first and the second detectors is constant 
and equal to 7r/2. Empty circles show the intensity distribu- 
tion (the counting rate R3 of detector 3). Solid line is the 
theoretical fit with Eq. (|4|. 



FIG. 6: Experimental setup for observing third-order polar- 
ization interference of classical light. GP, Glan prisms; S, slit; 
EOM, electro-optic modulator; BSl and BS2, beam-splitters; 
Oi-3, orientation of linear polarizers; Al-3, apertures; Dl-3, 
avalanche photodiodes; TCC, triple coincidence circuit; NDF, 
neutral density filters. 




FIG. 7: Polarization interference in the third order in the 
intensity observed by synchronously rotating two of the three 
Unear polarizers in opposite directions. Solid line shows the 
theoretical fit with Eq. (|16|) . Dashed line shows the fit taking 
into account a slight mismatch of the phases, 012 = 032 = 
7r/6, 013 — 
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FIG. 8: (a) Interference pattern recorded in a single laser 
pulse. The rectangle shows the area over which the analysis 
was performed: first, the image was averaged over the ver- 
tical coordinate (y) and then the resulting one-dimensional 
distributions were processed according to Eqs.([T7|. (b) The 
intensity distribution I(x) averaged over 500 shots. 
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FIG. 9: Third-order (a,c) and fourth-order (b,d) interference 
for coherent sources (a,b) and thermal sources (c,d). Dashed 
Une: theoretical fit with Eq. (|4| (a), Eq. (0 (b) and Eq. {[T3j 
(c). 



